How does torsional rigidity affect the wrapping transition of a semiflexible chain 

around a spherical core? 
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We investigated the effect of torsional rigidity of a semiflexible chain on the wrapping transition 
around a spherical core, as a model of nucleosome, the fundamental unit of chromatin. Through 
molecular dynamics simulation, we show that the torsional effect has a crucial effect on the chain 
wrapping around the core under the topological constraints. In particular, the torsional stress (i) 
induces the wrapping/unwrapping transition, and (ii) leads to a unique complex structure with 
an antagonistic wrapping direction which never appears without the topological constraints. We 
further examine the effect of the stretching stress for the nucleosome model, in relation to the unique 
characteristic effect of the torsional stress on the manner of wrapping. 

PACS numbers: 



I. INTRODUCTION 



The structure of DNA wrapping a cationic protein is 
called nucleosome, which is a basic unit of chromatin in 
living cells In physics, a spherical or cylindrical core 
and a semiflexible chain are adopted as a model to fea- 
ture nucleosome structure 043 ■ I n contrast to a flexible 
chain, a semiflexible chain wraps the core orderly [Io| . 
Under usual aqueous conditions, the persistence length 
l p ~ 50 nm of DNA is large compared to its molecu- 
lar thickness d ~ 2 nm. Therefore DNAs with contour 
length L>1 /im is characterized as a typical semiflexible 
chain: for example, the contour length of eukaryote DNA 
is longer than 500 /im. DNA wraps a histone core, whose 
diameter is 11 nm, 1.75 times orderly. To reveal struc- 
tures and mechanisms of nucleosome, some groups study 
by pulling nucleosome [Ir3 - [l5T ] . It was found that unwrap- 
ping is generated in a stepwise manner, i.e., the chain 
unwraps the core twice to once, then once to zero when 
the chain is pulled. These wrapping transitions are con- 
trolled by several factors, an interaction between them, 
the chain stiffness, and the sphere size. Recently, it is 
revealed that the other factor is also important to nucle- 
osome structures and mechanisms. Yanao et al. Q shows 
that DNA wraps protein core left-handed-like considering 
the coupling between bending and twisting rigidities owe 
to the chirality of the right-handed double helix struc- 
ture of DNA. In the system of a chain wrapping a core, 



topology is one of the important problems to reveal the 
mechanism and the structure of nucleosome. 

In prokaryote cell, DNA is a circular chain and topo- 
logically constrained. In these cases, topological num- 
ber Lk = Wr + Tw is conserved [l6[ . Wr is calculated 
through the chain conformation. Tw is total twisting 
of the chain. The effects of topological constraints and 
torsional rigidity on the conformation have been stud- 
ied in circular polymer [l7l - l20l |. In eukaryote cell, DNA 
is compactly stored by forming a hierarchical structure. 
In the lowest level, nucleosomes are packed into the 30 
nm chromatin fiber. In interphase cell, this fiber is fur- 
ther organized into Lampbrush chromosomes, which con- 
sists of highly condensed chromatin parts and extended 
chromatin loop parts emanating from condensed parts. 
Most DNA parts are generally not expressed, which are 
highly condensed. On the other hand, the chromatin 
loops are expressed actively [ll|. The chromatin loop 
is fixed both ends and topologically constrained. In eu- 
karyote, the structure of DNA is also influenced by the 
torsional rigidity. There are some researches to study 
how the topological constraints influence the conforma- 
tion of nucleosomes or chromatin fiber f2"2hj24j ] . However, 
there is an unsolved question how torsional stress under 
topological constraints affects the wrapping manner of 
the nucleosome. In this article, we will study the relation 
between torsional effect and the stability of the wrapping 
manner of nucleosome. We use a coarse-grained torsional 
model 17] to reveal the simple relation. We use a semi- 
flexible chain and a spherical core model inspired by nu- 
cleosome. 
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We report the nontrivial relation between torsional ef- 
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feet and the wrapping transition of semifiexible polymer 
based on a systematic Molecular Dynamics simulation. 
Our paper is organized as follows: In Sec. II, we intro- 
duce the model of a core and a semifexible chain which 
is adapted in the simulation. In Sec. IE, we show our 
simulational results. In Sec. IV, we discuss the wrapping 
transition in terms of free energy. Further results are in- 
terpreted in terms of free energy. Finally we conclude in 
Sec. V. 



II. METHODS 

To investigate the wrapping and unwrapping transition 
of semifiexible polymers, we carried out off-lattice Molec- 
ular Dynamics simulations in three-dimensional space. 
The model of one polymer and one core is essentially 
the same as that studied previously in Ref . 0, [l(| • Tor- 
sional potential is introduced following methods proposed 
in Ref. [Vi\ . The polymer is described as a beads-spring 
with positions r,; bond vectors = (r^+i — r,-)/|rj_|_i— r^|. 
The position of a spherical core is denoted by r c . To ac- 
count for the material twisting, two normal vectors and 
Vj are set at the gravity center of monomers rj and obey 
Vj = Ui x and • = 0. The potential energy of the 
system is represented by the following five terms; 
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set e histone = 5.0. The size of core a hlstone = 2.0a and 
CT/ = {&histone + a)/2. We set N = 60 polymer beads and 
one core. In this article, a, l p , and ahistone correspond 
to 5, 50, and 10 nm, which are good approximations to 
DNA and histone core (see introduction). 

Even though a tightly wrapped complex is formed, the 
core slides along the chain and prefers positioning at the 
chain end [J]. In order to eliminate possible end effects 
of the chain and to get clear-cut configuration, we fix the 
center of the core as like Ref. [l(| • 

The force f and the torque r is calculated from Eq. (1) 
- (5) and Eq. (4), respectively. The monomers obey 
the stochastic dynamics described by Langevin equation 
without momentum term 



71 dt~ fl +e ' 



dt 
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where is the Euler angle of pure rotation, r\ is the 
drag coefficient for drift, and r)R is the drag coefficient 
for rotation. The constant r = r\a 2 jk B T is chosen to 
be the unit for the time scale. We set the time step as 
dt = 2.5 x 10~ 5 t. The Brownian force & and torque £f 
satisfy the fluctuation dissipation theorem 



< &(*)&(*') >= 6k B Tr,6 l3 5(t - ti) 



(8) 



N 



Ubend _ kg .. 

k B T ~ 



< £ >= 2A B T^*(t - tr) 

(2) We set 77 = 1.0 and r] R = 0.213 

Let us introduce the following order parameter 
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where 9i is the angle between adjacent bond vectors, 0^ + 
7, is twist angle, which satisfies sin(ct; + 7;) = (vj • fj+i — 
fi • v i+ i)/(l + ^ • u i+ i) and cos(ai + 7;) = (f, • + 
Vj • Vi + i)/(l + U; ■ u,-(-i). The monomer size a and k B T 
are chosen as the unit length and energy, respectively. 
We set the spring constant kb on d — 500, the bending 
elasticity kg = 60, which corresponds to a persistence 
length l p ~ 10<7, and torsional energy k tor = 60. The 
excluded- volume effect is included in the Lennard- Jones 
potential Ulj, we set e = 0.20. For attractive energy and 
excluded- volume effect between a polymer and a core, we 
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where p(i) is an indicator of the pair contact: p(i) = 1 
if the number of monomers, which satisfies |r,; — r c | < r; 
and otherwise pij = 0. In the following discussion we set 
n = 2.5a. This quantity represents the degree of contact 
between the polymer and the core. As we shall see later, 
P is directly related to the wrapping number Nyy which 
measures how many times the chain wraps the core. 

In this paper, we study the chain conformation under 
three different manner of constraints; (i) twist one end 
(the other end is fixed), (ii) fixed both ends, and (iii) 
pull one end (the other end is fixed). In these handlings, 
the position of five monomers at the end (i"tv_4 to r^r) 
are fixed and '5 of those are also fixed. The positions of 
the other end monomers (ri to r$) are fixed in the case of 
(i) and (ii). In the case of (iii), those are fixed to y and z 
directions and pulled to x direction. $1 to ^ 5 are twisted 
in (i) and fixed in (ii) and (iii). We set = (60,0,0), 
rjv-i = (59,0,0), y N - 2 = (58,0,0), rjv_ 3 = (57,0,0), 
and rjv-4 = (56,0,0). The monomers of ri to are 
set at (y,z) — (0,0). The core is set at r c = (45,2,0). 
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We manipulate 5 monomers of both ends for the reason 
to prevent the chain crossing over the chain end and for 
convenience to twist and fix. 

We introduce topological number Lk, Wr, and 
Tw In this paper, both ends of the chain are fixed, 
Lk is conserved in the case of (ii) and (iii). In the case 
of (i), we decrease or increase Lk linearly. Lk is defined 
as follows 

Lk^Wr + Tw (II) 
Wr is defined as follows 

AvWr= ( [ (dri x dr 2 ) ■ r 1)2 /|r li2 | 3 (12) 
Jc Jc 

where ri and r 2 are the points passing along the closed 
curve C, ri 2 = r 2 — ri. In the calculation, we use 
Ref. 0. Tw is defined as follows 

N 

i=l 

We checked Lk every 100 steps in order to look out the 
unphysical topological breaking. 



III. RESULTS 

We first study how the chain wraps the core with twist- 
ing; Lk (control parameter) is decreased linearly from 
Lk = 2 to Lk = in 2 x 10 9 steps (t = 5 X 10 4 r) without 
topological breaking. The end to end distance is 33.5<r. 
Fig. [T] shows (A) evolution of topological numbers Lk, 
Wr, and Tw versus time and (B) snapshots of typical 
complex structures at t = 0, 2.5 x 10 4 , and 5 x 10 4 , re- 
spectively. The initial linking number Lk=2 indicates 
that the chain accommodates itself to the state with no 
twisting penalty by wrapping around the core twice in a 
left-handed fashion (Wr = 2) (Fig.Q](B) (a)). When the 
linking number is decreased, the complex first responds 
by the negative twisting while the writhing number kept 
almost fixed. This linear response regime is terminated 
at Lk ~ 1.5, around which the complex exhibits the 
global structural change, i.e., unwrapping transition from 
Wr ~ 2 to Wr ~ 1, and releases the torsional stress. Be- 
cause of the smallness of the system, the transition is not 
very sharp, but there are a finite range of the control pa- 
rameter (Lk), where the bimodal distribution is realized. 
The state with Wr ~ 1 and Tw ~ is stable around 
Lk = 1 (Fig. ED (B) (b)). By further decreasing Lk, the 
complex exhibits the second global structural transition 
at Lk ~ 0.5 to the state with Wr ~ 0. This, however, 
does not correspond to the unwrapping, but rather to 
the wrapping transition. Fig. [2] plots the order parame- 
ter P as a function of Lk, which clearly demonstrates the 
second structural transition as a wrapping. The inspec- 
tion of the snapshot indicates that each of two turns is 



characterized by the opposite handeness, and this antag- 
onistic wrapping results in Wr = 0. To realize such an 
antagonistic wrapping, there must be a "loop" in which 
the chain segment cannot be attached to the core (des- 
ignated by an arrow in Fig. Q] (B) (c)). We have con- 
firmed that its energetic stability is almost the same as 
that of the natural wrapped state; the total energy and 
the elastic bending energy of the chain in the case of 
Lk = are almost same as (little less than) the energy 
in the case of Lk = 2. The mean total internal energy 
u = U/(N - 10) (U is sumation of eq. (J) to (0), the 
mean elastic bending energy u bend = U ben d/(N - 10), 
and the mean adsorption energy u a d — U a d/(N — 10) 
are given as u = 1.264, u bm d — 1.066, u a d — —0.981 
in the case of Lk = and u = 1.269, u oen d — 1.086, 
Uad — —1.024 in the case of Lk = 2. It is interesting to 
note that although the presence of such a loop is an un- 
natural form, i.e., never observed in the system without 
topological constraint, the fluctuation of P is more sup- 
pressed in the antagonistic wrapping state (Nw = 2 and 
Wr — 0) than in the regular wrapping state (Nw = 2 
and Wr = 2). This structural stability of the complex is 
probably rendered by the steric hinderance of the chain 
segments which cannot cross one another (topological ef- 
fect of the chain). 

We study how the chain wraps the core with Lk. The 
both ends of the chain are fixed. The end to end distance 
is 33.5cr. Fig. [3] shows the distribution of wrapping num- 
ber versus Lk. In k tor = (without torsional effect), the 
peak is at Nw = 2, which indicates that the chain wraps 
the core twice. In Lk = 0, Lk = 0.25, and Lk = 0.5, the 
peak is at Nw = 2. These results are alomost same as 
the result of Fig. [3] In Lk — 0.75, there are two peaks 
at Nw = 1 and Nw > 2, which indicates that the chain 
wraps the core once and more than twice. In Lk ~ 0.75, 
the chain wraps and unwraps the core (see Fig. [2]) . In 
Lk — 1, the peak is at Nw — L The chain wrapping 
the core once is more favorable than wrapping the core 
twice although the chain loses the adsorption energy. In 
Lk = 1.25, there are no large peak; the probability is 
wide from Nw = 1 to Nw = 2. In Lk = 1.5, Lk = 1.75, 
and Lk = 2, the peak is at Nw — 2. These results are 
also same as the result of Fig. [2] 

These results indicate that Lk decides the stable wrap- 
ping number, how many times the chain wraps the core. 
In Lk = to Lk = 1, the peak shifts from wrapping the 
core twice to wrapping the core once. Then, in Lk = 1 
to Lk — 2, the peak shifts from wrapping the core once 
to wrapping the core twice. 



IV. DISCUSSIONS 

A. Internal energy 

A stiff circular DNA molecule is known to exhibit a 
buckling transition to a supercoiled state upon the in- 
crease in Lk [1^1 . Here the instability is caused by the 
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(A) 



Nw =1 Nw =2 




(B) 

(a) Lk=2, Wr=2, Tw=0 (t=0 ) 
Nw=2 




(b) Lk=1, Wr = 1, Tw=0 (t=2.5x10 ) 
Nw=1 / 

-9 



(c) Lk=0, Wr=0, Tw-0 (t=5x10 ) 
Nw=2 




FIG. 1: Quasi-static time evolution of the topological num- 
bers and the chain conformation. (A) Changes of Wr (red) 
and Tw (green) numbers accompanied by linear decrease of 
Lk = 2 to versus time steps. The data is an average of ten 
processes. (B) Typical snapshots of the chain and the core: 
(a) Lk = 2 (b) Lk = 1 (c) Lk = 0. 




FIG. 2: The wrapping number Nw versus Lk; Nw = 2 
(P=22.5) and Nw = 15 are described by solid and dotted 
lines respectively. The chain is twisted from Lk = 2 to lin- 
early with time (red line). The blue line is an average of ten 
processes. 



we now present a phenomenological theory to describe 
the properties of the core-chain complex with torsional 
effect. We refer the internal energy of Ref. [Tol. Il4j|. The 
internal energy of this case can be written as a sum of 
the adsorption, bending, and torsional energy terms. 



E(l) = E ad + E, 



bend 



E t , 



(14) 



where I is the chain length around the core. The energetic 
gain due to the chain adsorption E ad is written as E ad = 
el, where e is adsorption energy density (e < 0). The 
bending energy Ef, end is written as 



E bend = nl/2R 2 



--L)|sin,4| (15) 



where k is chain rigidity, R is radius of the core (ra- 
dius curvature), / is a extensional force at the end of 
the chain, and A is described by A = 2irl/4irR. The 
second term in eq. (|15p. which appears only in the case 
that both ends of the chain are fixed by pulling them, 
represents the penalty in the bending energy near the 
core (see Ref. [l3|). It has peaks at I = nR, 3ttR, which 
makes the halfhearted wrapped state with non-integer 
Nw unfavorable. The torsional energy E tor is written as 



balance between the twisting and bending energies. In 
our case of the twisting of the wrapped complex, the sim- 
ilar competition would be conceivable, but with an im- 
portant difference that the effective bending modulus can 
be regarded "negative", i.e., the chain favors the wrap- 
ping escorted by the adsorption energy, in the wrapping 
state (see discussion below). With such a point in mind, 



E,, 



W) 2 (^) 2 

A Lc 



2 Lc 



(16) 



where Lc is a contour length of the chain associated with 
the wrapping. Wr is related to the wrapping number 
Nw- We can write Tw — Lk — Wr ~ Lk — Nw in 
the regular wrapping state. Nw is roughly estimated as 
Nw — 1/2ttR. In the case of the antagonistic wrapping 
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(a) ktor=0 



Wr = \l- AttR\/2hR. Finally we get 





(b) Lk=0 to 1 



Lk=0 



(c) Lk=1 to 2 



Lk=1 



Lk=0.25 



Lk=1.25 



Lk=0.5 



Lk= 



Lk=1 




Lk=1.5 



= 1.75 



FIG. 3: Fhe distribution of the wrapping number with Lk: 
(a) Without torsional stress, (b) At Lk = 0, 0.25, 0.5, 0.75, 
and 1. (c) At Lk = 1, 1.25, 1.5, 1.75, and 2. 



state, 1 < N w < 2 and \Lk\ < 1, we use simulational 
results that the chain wraps the core twice but Wr = 0; 



(1< N w < 2, \Lk\ < 1) 
(17) 

We set e = -5.2, R = 1.0, k = 10, k tor = 10, / = 0.5, and 
Lq = 50. Fig. [4] shows the internal energy; (a) without 
torsional effect (b) with torsional effect as a function of 
Nw and (c) the value of the internal energy at Nw = 0, 
1, and 2 versus Lk. In this discussion, we define that the 
wrapping number Nw increases from to 2 linearly with 
a increase of I from to AirR. Without torsional effect 
(see Fig. 2](a)), the most stable is Nw = 2 and the stable 
state is N w = 1; E(l) Nw = 2 < E(l) Nw=1 < E(l) Nw=0 . 
In Lk = to Lk = 1, the internal energy at Nw = 2 
increases (see Fig. [U(b), (c)). In Lk = 1 to Lk = 2, the 
free energy at Nw = 2 decreases. The internal energy at 
Nw — 1 decreases with Lk = to Lk = 1 and increases 
with Lk = 1 to Lk = 2. These results indicate that the 
chain wrapping the core twice Nw = 2 is the most stable 
state in Lk ~ and 2 but not in Lk ~ 1. On the other 
hand, the chain wrapping the core once Nw = 1 is the 
most stable state in Lk ~ 1. Both terms E a d and Eb en d 
in eq. (14) have a linear dependence on kI/(2R 2 ) + el = 
e*l, so the wrapping transition takes place when e* < 

k/(2R 2 ) < \e\ in the absence of the torsional effect. 
Therefore, the torsional effect, if included, is a dominant 
factor in the total free energy of the system. The stable 
states are decided through decreasing torsional energy. 
We show a example here. In Lk ~ it is the most stable 
state to wrap the core twice. When the chain wraps the 
core twice, the adsorptional energy gain is high and the 
torsional energy is low because |Tiu| ~ (Lk = 0, Tw ~ 
0, and Wr ~ 0). When the chain wraps the core once, the 
adsorptional energy gain is low and the torsional energy is 
high because \Tw\ ~ 1 (Lk = 0, Tw±l, and Wr^l). In 
Lk = to Lk — 1, the stable state shifts from wrapping 
the core twice to wrapping the core once. Then, in Lk = 

1 to Lk = 2, the stable state shifts from wrapping the 
core once to wrapping the core twice. These discussions 
are consistent with the results of Fig. [3l 

We discuss the internal energy here because the en- 
tropy of this case should be almost zero: the free energy 
of this system is estimated as the internal energy. In 
opposite to this case, fluctuations and entropy are im- 
portant factors in the case of polynucleosome. 



B. force responses 

We study the stability of the wrapping around the core 
against the stretching. Fig. [5] (a) is tensional forces of 
pulled end (at r 5 ) versus L = |r 5 — rAr_ 4 | (like end-to- 
end distance). Fig. [5](b) shows the number of monomers 
around the core P versus L with and without torsional 
effect. In k tor = (without torsional effect), there are 
two peaks at L ~ 30, 40; it is a same result in Ref. fToj j . 
At L ~ 30, the chain unwraps the core twice to once; 
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(a) ktor=0 
^ 16 n 




- -4 I , , , J 

0.5 1 1.5 2 
Nw 



(b) 

Lk=0 — Lk=0.5 ~ Lk=1 — 
Lk=1.5 -- Lk=2 

16, , , , 



J2 12 




0.5 1 1.5 2 

Nw 



Nw =0 — Nw =1 — Nw =2 — 




0.5 1 1.5 2 
Lk 



FIG. 4: The internal energy of wrapping number Nw versus 
Lk by theoretical results, (a) in the case of ktor = (without 
torsional effect), (b) Lk — 0, Lk = 0.5, Lk = 1, Lk — 1.5, 
and Lk — 2 are drawn by black, green, blue, dashed red, and 
dashed yellow, respectively, (c) The value of internal energy 
at Nw = (black), Nw = 1 (red), and Nw = 2 (blue) versus 
Lk. 



Nw — 2 changes to 1. At L ~ 40, the chain unwraps the 
core once to zero; Nw changes to 0. In Lk = 0, there is a 
large extensional force at I ~ 30 because the chain wrap- 
ping the core twice is stable; E(1)n w= 2 < E(1)n w —%. On 
the other hand, there is a small peak at I ~ 37 because 
the potential barrier between E(1)n w= i and E(1)n w =o is 
small. The chain easily unwraps the core. In Lk = 0.5, 
there are two peaks at L ~ 30, 40. Nw is almost 
same behavior as the case of fc tor = 0; free energy are 
E{1)n w =2 < E(l) Nw= i < E(l) Nw = . In Lk = 1, there 
is a large extensional force at I ~ 37 because the chain 
wrapping the core once is stable. There is no peak at 
unwrapping the core twice to once because the chain 
unwraps the core easily: the potential barrier between 
E(1)n w =2 and -E(Z)jVw=i is small and the chain wrap- 
ping core twice is metastable (E(1)n w =2 > E(1)n w= \ < 
E(1)n w —o)- In Lk = 2, there is a large extensional 
force at L ~ 34 because the chain wrapping the core 
twice is stable. These is a very large extensional force 
at L ~ 42 because the chain not wrapping the core is 
much more unstable than the chain wrapping the core 
once; E(l) Nw=2 < E(l) Nw=1 < E(l) Nw = - In the case 
of \Lk\ > 0, E(1)n w= o is large, because Nw — implies 
Wr ~ in the present condition, thus the free energy 
quadratically increases with Lk; E(l) ~ E tor — (Tw) 2 ~ 
(Lk) 2 (see Fig.@](c)). It is noted that although the shape 
of real histone core is close to a cylinder, our spherical 
core model should capture the essential features in the 
wrapping-unwrapping transition. 

In this manipulation, we fixed five monomers of both 
ends in order to prevent topological breaking. This cor- 
responds to the case of adsorbing DNA end onto a tip 
of AFM, a micro bead, and the surface of a glass. In 
opposite to manipulate one monomer (to fix DNA at one 
point), it needs stronger force to pull monomers (to fix 
DNA at several points). There is a kink between a free 
monomer and a fixed monomer. It causes extra bending 
energy, which is estimated about 4/csT in this calcula- 
tion. 



V. CONCLUSIONS 

In this study, we have used a model inspired by nu- 
cleosome; a single semiflexible chain wraps a spherical 
core. The stable wrapping number, how many times the 
chain wraps the core, was studied under topological con- 
straints. With twisting once (quasi-static process), a in- 
crease in Lk (topological number) from to 1, the sta- 
ble states are shifted from the wrapping around the core 
twice to once. With twisting once more, a increase in 
Lk from 1 to 2, the stable states are shifted from the 
wrapping around the core once to twice. Internal en- 
ergy of the system including the torsional rigidity shows 
the same results as simulational results. The torsional 
energy is dominant in the case that the bending energy 
and the adsorption energy are erased each other; they lin- 
early depends on the length around the core. It restrain 
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(a) 



(b) 



ktor =0 


Lk-0 




Lk=0.5 


Lk=1 — 


- Lk-2 





20 25 30 35 40 45 25 30 35 40 20 25 30 35 40 45 
L L L 



Nw=0 Nw=1 



Nw =2 
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L 



25 30 35 40 
L 
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FIG. 5: Unwrapping process by stretching. The data are averages of ten processes, (a) Extensional force versus the distance 
of L (end to end distance); k tor = (without torsional effect), Lk = 0, Lk = 0.5, Lk = 1, and Lk = 2 are drawn by black, 
red, green, blue, and dashed red respectively, (b) Nw versus L. The wrapping numbers Nw = 2 (solid line), 1 (dotted line), 
(dashed line) are defined by P — 22.5, 15, 6 respectively. 



the wrapping conformation of the chain. The wrapping 
number, which is related to Wr, and the torsional stress, 
which is related to Tw, are coupled as Lk = Wr+Tw. In 
order to decrease torsional stress, the wrapping number 
needs change. The stability of the wrapping number is 
dependent on Lk through torsional stress. 

Finally, we mention the unwrapping process by stretch- 
ing. In the case increasing torsional stress after unwrap- 
ping, there is a large extensional force because the chain 
changes to less stable states. On the other hand, there is 
a small peak in the case of dereasing torsional stress after 
unwrapping because the chain easily unwraps the core; 
the chain is metastable state and the potential barrier is 
small. These results are easily verified by experiments. 
Recent experimental technic has developed and it is pos- 
sible to control torsion. For example, there is a study 
to twist and pull DNA 27]. We think that it is possi- 



ble to apply this experiment to pulling nucleosome with 
torsional constraints. We hope that our study will stim- 
ulate further experimental and theoretical development 
in the mechanical stability of chromatin and its link with 
biological functions. 
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